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Nonlinear localized magnetic excitations in one dimensional magnonic crystal is investi¬ 
gated under periodic magntic field. The governing Landau-Lifshitz equation is transformed into 
variable coefficient nonlinear Schrodinger equation(VCNLS) using sterographic projection. The 
VCNLS equation is in general nonintegrable, by using painleve analysis necessary conditions for the 
VCNLS equation to pass Weiss-Tab or-Came vale (WTC) Painleve test are obtained. A sufficient 
integrability condition is obtained by further exploring a transformation, which can map the VCNLS 
equation into the well-known standard nonlinear Schrodinger equation. The transformation built a 
systematic connection between the solution of the standard nonlinear Schrodinger equation and VC¬ 
NLS equation. The results shows the excitation of magnetization in the form of soliton has spatial 
period exists on the background of spin Bloch waves. Such solution exisits only certain constrain 
conditions on the coefficient of the VCNLS equation are satisfied. The analytical results suggest 
a way to control the dynamics of magnetization in the form of solitons by an appropriate spatial 
modulation of the nonlinearity coefficient in the governing VCNLS equation which is determined by 
the ferromagnetic materials which forms the magnonic crystal. 
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I. INTRODUCTION 


The study of nonlinear magnetic excitations in 
terms of solitary waves and solitons in ferromagnetic sys¬ 
tems have attracted much interest in the past several 
years[l-7]. The results reveal that the dynamics is gov¬ 
erned by Landau-Lifshitz equation which can be mapped 
to Nonlinear Schrodinger(NLS) family of equations[8]. In 
recent years, the studies on nonlinear systems with spa¬ 
tial periodicity has become a great topic of interest [9]. 
BEC in optical lattices[10,11], solitons in Photonic lat¬ 
tices [12] and periodic magnetic systems[13-15] etc., are 
the typical models among them. Motivated by these 
considerations, in the present paper we investigate the 
nature of excitation of magneization in one dimensional 
magnonic crysatal. Magnonic crysatal is a medium with 
spatially periodic variation of their magnetic properties 
in a definite direction. In the linear regime, observation of 
frequency band gap is well studied problem. The funda¬ 
mental feature of periodic magnetic structures is energy 
band gap in their spectrum of spin waves. The band 
gap represents a range of energy values in which spin- 
wave excitations are forbidden from propagating. The 
theoritical and experimental studies related to magnonic 
crystal mostly is devoted to linear phenomena. The in¬ 
vestigation of propagation of soliton in magnonic crystal 
are insufficent, only few specific studies in the field that 


shows the experimental and numerical simulation results 
based on one dimensional NLS equation. The bright and 
dark solitons were observed in yittrium iron garnet films 
with artificial periodicity[13-15]. Morozove et.al[16], in¬ 
vestigated the features of formation of the soliton that 
are similar to bragg solitons in the ferromagnetic one 
dimensional periodic structure using coupled mode the¬ 
ory. He et.al[17], studied the modulation instability and 
gap solitons in ferromagnetic films under periodic mag¬ 
netic field using multiscale expansion method. The ear¬ 
lier studies are based on a homogeneous ferromagnetic 
films and achieve periodicity by varying the thickness 
of the films or by applying spatially periodic magnetic 
field[16,17]. Here in this present study, we consider an 
infinite one dimensional magnonic crystal formed by pe¬ 
riodic array of distinct elements and study the impact 
of material parameter variation on the localized exci¬ 
tation of magnetization under periodic magnetic field.. 
The paper is organized as follows. In Sec.H we consider 
the one dimensional magnonic crystal model under peri¬ 
odic magnetic field and derive the dynamical equation. 
In Sec.HI The governing VCNLS equation is analyzed 
through painleve analysis to obtain integrability condi¬ 
tions and it is mapped into standard NLS equation using 
suitable transformation. In Sec.IV The results are pre¬ 
sented. 


2 


II. MODEL AND DYNAMICAL EQUATION 

We consider an infinite one dimensional 
magnonic crystal represented by a system of alternating 
uniform ferromagnetic layers of two different materials 
A and B as shown in Fig.(l). The layers have different 
values of exchange length parameter, and saturation 
magnetization,throughtout the sample. Let Jex,A 
and Jex,B be the exchange length for the ferromagnetic 
materials A and B respectively. Let Ms ,a and Ms,b be 
the saturation magnetization for the ferromagnetic mate¬ 
rials A and B respectively. Here z-axis is chosen normal 
to the plane of the layers. The equation of motion of the 
magnetization in the ID-magnonic crystal is governed by 
the following Landau-Lifshitz(LL) equation[18]. 

dMir,t) ^ x Aff{r,t), (1) 

where 7 is the gyromagnetic ratio and Hef / denotes the 
effective field. In general, the effective field Heff is the 
sum of several components includes the applied field, 
anisotropy field, demagnetization field and exchange field 
which are all dependent on space. 

Heff — Hq -|- Hani Hqx- ( 2 ) 

The first component Hq is the applied magnetic field 
which is inhomogeneous in space and applied along z- 
direction. The next component is the anistopy field, Hani 
which is given by 

Hani = Pix)M^Z. (3) 

The another component Hd arises entirely from the de¬ 
magnetizing field that corresponds to shape anisotropy. 

Hd{r, t) = X{NxMxX + NyMyi) + N^M^z), (4) 

In the case of an xy-film, the demagnetization field is 
given by Hd = XMzZ. 

Where A = — 1 . 

The exchange field Hex given by [19], 

Hex = Jex{x)V‘^M{f,t)), (5) 

where Jex{x) = ^ exchange length, A{x) 

is the exchange constant and Ms{x) is the saturation 
magnetization. Thus, the total effective field Heff takes 
the form, 

Heff = Ho{x)z -\- j3{x)MzZ + XM^Z + Jex{x)'V^M{f,t))- 

(6) 


Upon using the above expression for the effective field in 
Eq.(l), we get 

= -jM X [Ho(x)z + P(x)M^z + 
ot 

+Jexix)V^Mir,t)], (7) 

Ml + Ml + Ml = Me{xf. 



FIG. 1. Schematic drawing of one dimensional magnonic crys¬ 
tal and its coordinate system. A and B are two ferromagnetic 
materials. 


Jex{x), Ms{x) and P{x) are periodic functions with pe¬ 
riod equal to the magnonic crystal lattice constant a. 
The exchange length, Jex{x) is represented as. 


Jexi^^ “1“ — JexifHj — 


Jex,A 0 ^ X <C 


Jex,B n/2 X Qj. 

The saturation magnetization, Ms{x) is represented as, 

i\/f ( \ \ i\/T ( \ j^s,A 0 < X < n/2, 

Ms{x + a) = Ms{x) = < / 

I Ms,b cl/z < X < a. 

Similarly, the anisotropy constant, /3{x) is represented as, 

Pa 0 < X < a/2, 


P{x -h a) = P{x) = 


Pb < X < a. 


Since the Landau-Lifshitz equation is a continous 
equation, the material parameters should be represented 
in continous form. To make its continous which are 
expanded in a Fourier series in the form. 


T / \ ^ .‘2fl7TX. ^—> . .‘2t17TX. . . 

Jex[x) = Co + > CnCOS[ -) + > dnSin[ -). (8) 

a a 

n=l n=l 


After evaluating the coefficients of the Fourier series, we 
get the function Jex (x) in the form. 


JexifP) — ( 


Jex,A Jex,B 


n=l 


. 2n7TX . 
xsin{ -), 


(9) 


Jexix) = Jav + Y( —)(1 - (-l)"))sm(^), (10) 

nir a 

n=l 

The saturation magnetization Ms{x) is represented as, 

Ms{x) = ^ Ms,a + Ms,b -^ ^ 5](^)(1 - (-1)")) 

A ^^ niT 


n=l 


2n7rx 

xsini -), 

a 


( 11 ) 
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Ms{x) = Ms^av + V(-^)(1 - (-l)^))5zn(-). 

niT a 

n=l 

( 12 ) 

Similarly, the anisotropy constant 0{x) is represented as, 

I3a + I3b. . ^,A/3 , 2mTx 

P{^) = (-^-) + Z ^(—)(1 - (- 1 ) —)> 

Zi ihn CL 

n=l 

(13) 

m = Hav + E(—)(1 - (-l)"))sm(^). (14) 

^ nTT a 

n=l 


where, — Jex,A Jex^Bt 

A/3 = Pa — Pb and Jav, Ms^av and Pav represents the 
average exchange length, saturation magnetization and 
anisotropy constant value of two ferromagnetic materials 
A and B respectively. The material parameters varies 
gradually at the interface between the two ferromagnetic 
materials and the average values of the material param¬ 
eters represents, a value of material parameter of the pe¬ 
riodic ferromagnetic system at the exact centre interface 
point between the two ferromagnetic materials. 

The LL equation is a vector nonlinear partial differential 
equation, it is difficult to solve in its original form. By 
using sterographic projection, we transform the LL equa¬ 
tion into a nonlinear equation of a complex function. In 
the component form Eq.(7) becomes. 


^ = -^{Ho{x)My + p{x)M,My + XM,My 

+ JeAxKMyV^M,-M,V^My)}, (15) 

A = -j{-Ho{x)M, - I3{X)M,M, - XM,M, 

- Je,(a:)(MW'M, - M,W^M,)}, (16) 


dM, 

dt 

Defining, 


-AJe:o{x){M^V^My - MyV^M^)}. (17) 




Mx + iMy 

Msix) 


(18) 


^*{x t) = — _^ 

^ ^ ’ Ms{x) ’ 

where -0 is a complex variable and then we have, 

mz{x,t) = (1 - 


(19) 


( 20 ) 


Considering small deviations of magnetization from the 
equilibrium direction corresponding to << 1 and un¬ 
der the long wavelength approximation by keeping only 
the nonlinear terms of magnitude |'0p'0[2O] we obtain. 


-'‘-'>W'(S?))^ 


Ms{x) \ 
Ms^av ) 


Ho{x) 

Ms,av 


tp\ V’ 
( 21 ) 


Here the temporal, and spatial coordinates are rescaled 
by to = l/i'jMs^av) and lo = ^/iJav) respectively. 

Let 

/(*) = (^) (^) > 9ix) = (1 - m) (^) 

Then the Eq.(21) becomes, 

- f(x)^+ ^g{x)\ip\^'ip+ h{x)ip = 0 ( 22 ) 

The Eq.(22) is the Nonlinear Scrodinger equation with 
variable coefficients. 

When f{x) = g{x) = h{x) = constant Eq.(22) reduces 
to completely integrable nonlinear schrodinger equation 
which admits N-soliton solutions [21]. 

In the absence of cubic term, the above Eq.(22) is 
a linear periodic system and admits Bloch wave so¬ 
lutions. As mentioned earlier, in the linear studies 
observation frequency band gap is well studied problem. 
The bandgap represents the range of energy values 
in which spin wave excitations are forbidden from 
propagating. Erom the fabrication point of view, exper¬ 
imental studies were performed on periodic structures 
composed of only one constituent magnetic material 
and achieve one dimensional periodicity in homogeneous 
ferromagnetic films which includes micron size shallow 
grooves etched on yttrium iron garnet(YIG) films[22], 
one dimensional array of micron size metal stripes on 
YIG films [23] and by applying a periodic magnetic field 
of spatially varying strength[24]. Wang et.ah,[25,26] 
have achievd first bicomponent magnonic crystal experi¬ 
mentally which exhibit well defined frequency bandgaps. 
The experimental results revealed that the tunability 
of the magnonic band gaps can be achieved by varying 
the width of the component stripes or by varying the 
materials. 

In the presence of cubic term. It is completely nonlinear 
problem. The above VGNLS equation is in general 
nonintegrable, to solve the above equation analytically, 
in the next section integrability conditions are obtained 
by using painleve analysis. The integrability conditions 
are expressed in terms cofeiheients of the VGNLS 
equation. The above equation admits soliton solutions 
only when the integrability conditions are satisfied. 


III. PAINLEVE ANALYSIS AND 
INTEGRABILITY CONDITIONS 

Several tools such as Painleve analysis[27], Lax 
pair [28] and similarity transformation techniques [29] are 
available to solve VGNLS equation to obtain analytical 
solutions. Our analysis is based on the Painleve test 
for partial differential equations i.e., the Weiss-Tabor- 
Garnevale (WTG) test, which has been found to be a 
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successful tool for investigate the integrability of partial 
differential equations. In this section, we use WTC test 
to obtain an integrability condition for the VCNLS equa¬ 
tion and then under this condition, we look for a transfor¬ 
mation which converts the Eq. (22) to the standard NLS 
equation. In order to perform conveniently, we rewrite 
the Eq. (22) and its complex conjugate by replacing -0 
by a and -0* by b obtain, 

~ ^ - 0 (23a) 


M*) = (iSS) - ( 


Ms(x) 

^S, a V 


) 


Substuting this in eq.(37) and eq.(38) we get the 
integrability conditions as. 


and 


Jex (^) 
Jav 


Mlix)) 


(28) 


“*1? “ ^ + Kx)b = 0 

Where /(x), g{x) and h{x) are real functions. 


(23b) 


Ho{x) \ 
Ms,av ) 


i.e., 


(Ms{x)\ _ (sM'sHx) + 2Ms{x)M''{x) 

V Ms,av ) V 

(29) 


The next step is to seek solution in the form of 
Laurent series. 



a{x,t) = T,^Qaj{x,t)(j)°‘'^^ {x,t) 

(24a) 


b{x,t) = W=obj{x,t)r+^ix,t) 

(24b) 

Where ao, 6o 7^ 0 and a^, bj and (l){x,t) are analytic 
functions, a and jS are negative integers to be determined 
from the leading order analysis. 

Atr j=4, we obtain integrability conditions. 


+ 

to 

II 

0 

(25) 

i.e.. 


(26) 


'■<^'71) [f-lf] 

(27) 


Hojx) _ f Ms{x) \ 2 ( -SMsHx)+2Ms{x)M''{x) 

Ms,av \Ms,avJ ( M|((r) 

(30) 

The above integrability conditions are consistent with 
Ref.[10] and Ref.[27]. The exchange constant, Jex{x) 
which is inhomogeneous in space related to the satura¬ 
tion magnetization Ms{x) is given by Eq.(39) and it is 
first integrability conditions for Eq. (22) to be integrable. 
Erom Eq. (41) it is noted that form of the periodic ap¬ 
plied magnetic field is determined by the ferromagnetic 
materials which forms the magnonic crystal. The ap¬ 
plied magnetic field is periodic in space and its period¬ 
icity forms a periodic potential for the spin waves. The 
form of the potential is related to saturation magnetiza¬ 
tion Ms{x) which is given by Eq.(41) and it is second 
integrability conditions for Eq. (22) to be integrable. 
Eurther, we have to construct the solution of VCNLS 
equation by using a transformation which converts the 
Eq. (22) into a standard NLS equation. We look for the 
transformation of the form[10], 

'Ip{x,t) =r{x)q{X{x),T{t)) (31) 


Where k is an integration constant. By employing, 
Painleve method for the governing VCNLS Eq.(22) and 
we obtain integrability conditions. 

Here, 

/(-) = (^) (SiS). = (‘ -(ss) 

The difference in anisotropy constant values [3 between 
the constituent materials tends to slight and the effect of 
this inhomogenity is minor. Hence assume this compo¬ 
nent to be negligible. 

Then, 


where X = X{x) and T = T(t) and r{x) are the real 
functions to be determined. 

Substuting Eq. (42) into Eq. (22), We get set of following 
equations 


f{x)r:ox + h{x)r(x) = 0. 

(32) 

2vxX^ 4“ xix^X^^ — 0. 

(33) 

Tt = f{x)Xl = g{x)r{xf 

(34) 


and also, 

iqr - qxx + -^\Qf Q = 0 - 


fi’i = (^) (SiS)-= («S) 


(35) 



5 


By using the integrability conditions, the above equations 
are solved and gives 


r{x) = 



(36) 


X{x) = J g{x)dx and T(t) = t (37) 

Then, we obtain solution of Eq. (22) by known solution 
of standard NLS equation, q{X,T). 


= y -^q{X(x),T{t)) (38) 

Here, q(X,T) is the solution of standard NLS equation 
Eq.(46) and there exists several methods to solve the 
standard NLS equation such as classical 1ST, DBT, Hi- 
rota bilinear method etc.,. In this paper, Hirota bilinear 
method is used to construct the dark one soliton solution. 

q{X{x),T{t)) = TfCl - 2iC2ta.nh{C2{X - CIT))] 

v2 

X exp + 4C2^)T] (39) 


V. RESULTS 


Erom Eq.(49), we obtain 

= -j^=[Cl - 2zC2tanh(C2(X - CIT))] 

V2^(x) 

X exp (Cl^ + 4C2^)T] (40) 

where the parameters Cl and C2 corresponds to the 
velocity and depth of the dark soliton. 


g{x) = 1 + T] 


n=l 


Ms,A — Ms,b \ 
Ms,avmT J 


(1 - (-l)”))sm( 


2n7rx 

a 


(41) 

where Ms^a and Ms^b are the saturation magnetization 
for the ferromagnetic materials A and B respectively. 

Ms,av = ^ is the average saturation mag¬ 

netization value of two ferromagnetic materials A and B. 
The material parameter Ms{x) varies smoothly at the 
interface between the two ferromagnetic materials and 
the average value represents the saturation magnetiza¬ 
tion value of the periodic ferromagnetic system at the 
exact centre interface point between the two ferromag¬ 
netic materials. 

As we mentioned earlier, Landau-Lifishitz equation is a 
continous equation which describes the equation of mo¬ 
tion of the magnetization in a ferromagnetic medium. 
Here we consider a periodic ferromagnetic system in 
which the material parameter varies periodically. In or¬ 
der to make it continous and incorporate into LL equa¬ 
tion here we use fourier series to represent the periodic 


material variation into continous form. In the above 
equation, g{x) represents continous form for the varia¬ 
tion of saturation magnetization value at each points in 
the periodic ferromagnetic system having equal widths. 
By using different ferromagnetic materials to form a pe¬ 
riodic ferromagnetic structure, form of the g{x) changes 
accordingly. 

Erom -0 we obtain the components of magnetization, 

+ (42) 

my{x,t) =( ^ (43) 

= (44) 

The results indicates that the amplitude of the soliton 
solution depends on the nonlinearity coefficient g{x)^ 
which means that the soliton can be spatially modulated 
and which admits several interesting spatial phenomena. 

Case (i): 

When Ms,a=Ms,b^ then g{x)=l it represents homo¬ 
geneous ferromagneic system which is governed by 

standard nonlinear Schrodinger equation and admits 
soliton solution which propagates in a homogeneous 
background shown in Eig. (2). 


\x/\^ 




FIG. 2. Soliton propagates in a homogeneous background 

Case (ii): 

When Ms,aJ^J^s,b^ then 


9{x) 


1 + E 

n=l 


/ Mg,A — Ms,b \ 

V Ms,avmT J 


(1- (-l)”))sm( 


2n7TX 

a 


)• 


(45) 

g{x) = 1 + Insind^^^^). (46) 

a 

n=l 

Where In is the control parameter which determines the 
nature of the magnonic crystal. 
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Case (a): 

Consider a magnonic crystal system formed by periodic 
array of distinct ferromagnetic elements of iron,(Fe) and 
cobalt,(Co). The magnonic crystal is formed with lattice 
constant, a is 500nm by choosing the width of the each 
layer as 250nm. The spatial magnon density profile, 
|Tp and ruz components in 3D and 2D of magnonic 
crystal for the combination of 250Fe/250Co are shown in 
Fig.(3a-3d). Its corresponding periodic applied magnetic 
field which is the condition for integrability from Eq.(41) 
is shown in Fig. (4). The excitation of magnetization in 
the form of spatially periodic localized modes is exists 
in the oscillatory background with structure similar to 
the form of spin Bloch waves [30]. 

Lattice Constant of magnonic crystal a=500nm 
25QFe/25QCo _ 


Material 

Saturation 
magnet iza- 
tion Ms (10® 
A/m) 

Exchange length, Jex (nm) 

Ms, A, Fe 

1.752 

3.30 

Ms,b, Co 

1.445 

4.78 


Ms,av= 1.598x10^ A/m and Jex,av= 4.04nm 

Case (b): 



FIG. 3. (a) and (b) The spatial profile of magnon density |T|^ 
in the form of soliton on the background of spin Bloch waves 
in 3D and 2D. (c) and (d) rUz component in 3D and 2D for 
250Fe/250Co. 


and rUz components in 3D and 2D of magnonic crystal 
for the combination of 250Co/250Py shown in Fig. (5) 
and its corresponding periodic applied magnetic field 
which is the condition for integr ability from Eq.(41) 
are shown in Eig.(6). As in the previous case, the 
excitation of magnetization in the form of soliton has 
spatial period exists on the background of spin Bloch 
waves. The parameter g{x) determines the nature of 
the magnetic crystal and periodic applied magnetic field 
which act as periodic potential for spin waves to satisfy 
the integrability condition. The stability nature of the 
soliton solution is depends on these parameters which 
are material dependent and its completely under our 
control. By choosing the ferromagnetic materials which 
forms the magnonic crystal of our interest, it is possible 
to tune the spatially modulated amplitude of the soliton. 
This spatially modulated amplitude soliton solutions 
with oscillatory background describes the nonlinear 
localized exicitation of magnetization in one dimensional 
magnonic crystal. 

Lattice Constant of magnonic crystal a=500nm 
25QCo/25QPy _ 


Material 

Saturation 
magnetiza¬ 
tion Ms (10® 
A/m) 

Exchange length, Jex (nm) 

Ms, A, Co 

1.445 

4.78 

Ms,b, Py 

0.860 

7.64 


Ms,av= 1.152x10^ A/m and Jex,av= 6.21nm 





250 Fe/250 Co 
1.4 


^ 1.2 



0.6 


0 200 400 600 800 1000 

x(nm) 


FIG. 5. (a) and (b) The spatial profile of magnon density |T|^ 
in the form of soliton on the background of spin Bloch waves 
in 3D and 2D. (c) and (d) rUz component in 3D and 2D for 
250Go/250Py. 


FIG. 4. Periodic applied magnetic field which is the condi¬ 
tion for integrability from Eq. (41) for the combination of 
250Fe/250Go. 

Consider a magnonic crystal system made of cobalt,(Co) 
and permalloy,(Py) with lattice constant, a=500nm by 
choosing each width as 250nm. The density profile, |Tp 
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250 Co/250 Py 



0 200 400 600 800 1000 

x(nm) 


FIG. 6. Periodic applied magnetic field which is the condi¬ 
tion for integrability from Eq. (41) for the combination of 
250Co/250Py. 


V. CONCLUSION 

In this paper, by transforming the governing Landau- 
Lifshitz equation of one dimensional magnonic crystal 
into Variable coefficient Nonlinear Schrodinger equa- 
tion(VCNLS), we have investigated the dynamics of 
magnetization in one dimensional magnonic crystal 
and with aid of painleve analysis we constructed the 
soliton solution exists on the oscillatory background with 
structure similar to the form of spin Bloch waves. Such 
solutions exists in certain constraint conditions on the 
coefficients of the VCNLS equation. The results shows 
that the amplitude of the soliton solution has spatial 
period on the background of spin Bloch waves. The spa¬ 
tial distribution of the soliton profile determined by the 
free parameter g{x)^ which is depend on the saturation 
magnetization, Ms values of the ferromagnetic materials 
which are used to form the magnonic crystal. From the 
case studies it is observed that by varying the parameter 
g{x) nonlinear coefficient of the VCNLS equation for 
diiferent combination of the magnonic crystal, the form 
of the periodic applied magnetic field can be changed 
to satisfy the integrability condition and which is act as 
the potential barrier for the spin waves and accordingly 
the desirable amplitude modulation of the soliton can 
be achieved. 
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